Divide and conquer: resonance induced by competitive interactions 
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We study an Ising model in a network with disorder induced by the presence of both attractive 
and repulsive links. This system is subjected to a subthreshold signal, and the goal is to see how the 
response is enhanced for a given fraction of repulsive links. This can model a network of spin-like 
neurons with excitatory and inhibitory couplings. By means of numerical simulations and analytical 
calculations we find that there is an optimal probability, such that the coherent response is maximal. 



INTRODUCTION 

In nonlinear systems, the right amount of noise can 
amplify the response to a weak periodic signal. This con- 
structive effect of noise, known as stochastic resonance, 
was first proposed to explain the occurrence of ice ages 
[J, and has since found applications in such diverse 
areas as lasers Jsl, SQUIDS [J], or neurons 5], just to 
mention a few [a, lZ| ■ 

The mechanism of stochastic resonance involves a 
matching between the frequency of the external signal 
and a stochastic frequency induced by the noise. The pro- 
totypical system is one in a bistable potential subjected 
to a periodic modulation signal and to noise. A weak, 
subthreshold, signal cannot, by itself, make the system 
switch between its equilibrium states. When driven only 
by noise, the system is able to jump between the two 
minima, with a mean frequency which depends on the 
intensity of the noise as given by the celebrated Kramers' 
formula [sj . It is then possible to tune the noise intensity 
in order to match the noise induced jumping frequency 
with the frequency of the forcing. At this point, we have 
an amplified coherent response. Based on those simpli- 
fied ingredients, stochastic resonance has been applied to 
many areas and situations, and extended from systems 
with a few degrees of freedom to systems composed by 
many units. 

It was shown in reference [2] that diversity or hetero- 
geneity, in the form of quenched disorder, can play the 
same constructive role as noise. The authors considered 
a system composed by many coupled bistable units, sub- 
jected to an external weak periodic signal. Diversity is 
introduced as variability of a parameter that controls the 
relative stability of each bistable state of the potential. 
When the units are identical (and both states are equally 
stable for all units) the signal is subthreshold and, be- 
cause of the coupling, all units remain in the same state. 
As diversity increases, the signal becomes, for half of its 
period, supra-threshold for some of the units and forces 
those units to jump from their less stable state to the 
other. In the other half of the period, the signal becomes 
supra-threshold for a different set of units. The units 



which follow the signal pull the other units, to whom 
they are attractively coupled, and the collective effect is 
that a significant fraction of the units is able to respond 
to the external forcing. For this mechanism to be effec- 
tive, the units cannot be too diverse, because then some 
of them would offer too much resistance to follow the 
signal. But they can not be too similar either, because 
then there wouldn't a big enough fraction that can follow 
the signal. The optimal value of diversity is the one that 
makes the system more sensitive to the external signal. 
The collective effect can be understood as the result of 
the degradation of entrainment induced by diversity, a 
similar effect to that induced by noise This degra- 
dation results in the lowering of the effective potential 
barrier separating the two stable states of the collective 
system. The barrier can then be more easily overcome 
by the external forcing. One of the main points made 
in [9| is that it is the loss of entrainment that drives the 
resonance effect. This loss of entrainment can be induced 
by diversity [91, by noise (in the case of extended stochas- 
tic resonance (ill fl^). or by some other source. Along 
these lines, the role of the heterogeneous complex net- 
work topology in the amplification of external signals has 
been addressed in [l^l, and Chen et al. 14 1 have shown 
how structural diversity enhances the cellular ability to 
detect extracellular weak signals. The interplay between 
noise and diversity in an ensemble of coupled bistable 
FitzHugh-Nagumo elements subject to weak signal has 



been considered in [15 1 



The presence of both repulsive and attractive interac- 
tions is not unusual in systems with many units. The 
existence of inhibitory and excitatory connections in the 
brain neurons, or a society with friends and enemies are 
examples of such systems. The emergence of a coher- 
ent behavior in the absence of forcin g an d in the pres- 
ence of repulsive links was treated in [16l |. There it was 
shown that one can obtain a more coherent behavior, 
in the form of synchronized pulsing, by adding an opti- 
mal amount of long-range repulsive couplings in a mix- 
ture of excitable and oscillatory units described by the 
Hodgkin-Huxley model. In the same reference, a sim- 
ilar improvement of the internal coherence in an Ising 
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model with a simple majority-like dynamics in the pres- 
ence of long-range repulsive links was also shown. Also 
in [l7| . an intermediate amount of repulsive links was 
found to trig ger collective firing in an ensemble of active- 
rotators [l8| in the excitable regime. The role of diversity 
in heterogeneous excitable media was considered in [19| 
where the author demonstrates that diversity in a param- 
eter can cause the emergence of global oscillations from 
individually quiescent elements in a system of van der 
Pol-FitzHugh-Nagumo elements. The combined effects 
of noise and variability in th e sy nchronization of neural 
elements has been studied in [20|, while reference 21 1 un- 
veils the general mechanism for collective synchronized 
firing in excitable systems arising from degradation of 
entrainment originated either by noise, diversity or other 
causes. 

In this work we study a periodically forced system 
where the only source of disorder is competitive inter- 
actions and show that competition in the sign of inter- 
actions may also lead to a resonance effect. This reso- 
nance can be interpreted as an optimal transmission of 
the information carried by the external signal, in a kind 
of "divide and conquer" effect. The paper is organized as 
follows: In section we introduce the model and present 
the results of numerical simulations that show the ex- 
istence of the resonance effect; in section a mean-field 
approximation, which is able to reproduce qualitatively 
the simulations results is detailed; in section we discuss 
in detail the mechanisms that may lead to the resonance, 
both from microscopic and macroscopic points of view; 
finally, in section, we end with some brief conclusions 
and outlooks. 



MODEL AND RESULTS 
Model 

We consider a set of N spin-like (Ising) dynamical vari- 
ables fJ.i{t) which, at a given time can adopt one of two 
possible values, fXi = ±1. We will sometimes use the lan- 
guage of a magnetic system, but our aim is quite general 
and these states can represent, for instance, two differ- 
ent opinions (in favor/against) about a topic, the state of 
a neuron (firing/not firing), or several other interpreta- 
tions 



22l | . The variables are located on the nodes of a 
given network whose links represent interactions. We as- 
sign a weight coij to the link connecting nodes i and j and 



consider only the symmetric case u>i 



(or an undi- 



rected network). According to the discussion above, we 
let the weights take positive or negative values: — 1 
or LOij = —K with K > 0. The neighborhood of node 
i is the set V{i) of nodes j for which a connecting link 
between nodes i and j exists. 

The spin variables evolve according to the following 
dynamical rule: At time t one of the variables, say /x^, is 



chosen at random. The value of this variable is updated 
according to: 



sign 



E 



w.p. 



1 



sign [sin(rii))] 



|asin(f7i)|, 

w.p. |asin(rit)|, 

(1) 

(w.p. stands for "with probability"). In both cases, if 
the expression within square brackets is equal to zero, 
the variable does not change: /ii(t + r) = /ii(t). The first 
case represents a weighted "majority-rule" in which the 
state of the spin is determined by the sign of its local field 



h^{t) = 



fj-jit). The second case represents the 



effect of an external forcing of frequency f2 - the intensity 
a < 1 determines the rate at which the signal influences 
the dynamics of the variable /i^. The choice of the time 
step r — 1/N defines the unit of time as N updates. 
We consider both regular lattices (with k neighbors) and 
random networks of the small-world type. The latter 
are constructed in the usual way (23l |. Denoting by q 
the rewiring probability (percentage of short-cuts), the 
limit q — 1 corresponds to a random Erdos/Renye-type 
network, g = is a regular ring-network and intermediate 
values of q define a small-world network. We have also 
considered a square lattice in which a node is linked to 
the k = 8 nodes of its Moore neighborhood. In each 
case, links are assigned a strength with probability 
p or a strength 1 with probability 1 — p. In the case 
of a random network, the number of links (degree) ki 
of node i is a random variable with probability Pfc. and 
average {ki) — k. Denoting by k^ and k~ respectively 
the number of positive and negative links of node i, its 
degree is ki = kf -\- k~ and (k^) = (1 — p)k, {k~) = pk. 

It is worth noticing that, from the formal point of view, 
the majority-rule is equivalent to a heat-bath stochas- 
tic dynamics in the limit of zero temperature [i^]. The 



Hamiltonian is H ■ 



fiifij (the sum runs over 



all pairs of neighbors) and the majority-rule always leads 
to a configuration with less or equal energy. If all the 
weights LOij SLie positive, the ground states are /i^ — +1 
or /ii ~ —1, Vi, and these ground states are reached in- 
dependently of the initial condition. If there is a fraction 
of negative links, the system is of the spin-glass fam- 
ily. The (in general unknown) ground state can have 
many metastable configurations nearby and the use of 
the majority-rule may trap the system in one of them. 

As a way of quantifying the coherence of the global 
response to the forcing, we chose the spectral amplifica- 
tion factor R, defined as the ratio of the output to input 
power at the corresponding driving frequency 25| : 



R 



|((m(t)e-"*)) 



(2) 



where 



is a time average, ni(t) is the global response 
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FIG. 1: Spectral amplification factor R versus probability 
of repulsive links p. Parameters are; a = 0.15, S7 = 
K = 1. In the main graph, A'^ = 100 and symbols correspond 
to topologies: ring with = 10 neighbors (o), square lattice 
with fc = 8 neighbors in the Moore neighborhood (□), and 
random networks with average number of neighbors fc = 10 
and rewiring probability q = 0.2 (*) and q = 1 (a). In the 
inset, we chose the random network with g = 1, fc = 10, and 
different curves correspond to sizes N = 100 (a), 500 (<), and 
1000 (V). 

(system's magnetization): 

1 ^ 

and (...) is an ensemble average over network realizations, 
initial conditions and realizations of the dynamics. Large 
values for R indicate that the global variable m{t) follows 
the external forcing, while small values of R indicate a 
small influence of the forcing on the global variable. 



Simulation results 

The main result of this paper is that there is a reso- 
nance effect, a maximum of the amplification factor R, 
at an intermediate value of the probability of repulsive 
links p, as shown in figure [TJ The existence of this max- 
imum is also visible when looking at the amplitude of 
the oscillations of the global variable m(i) - figure [H For 
small p, m{t) oscillates with a small amplitude (of order 
a) around a value close to either -1-1 or —1. As p in- 
creases, one clearly notices that the amplitude increases 
dramatically and m.{t) oscillates around 0. As p increases 
even further, the amplitude of the oscillations decreases 
but the global variable still oscillates around 0. This res- 
onance effect appears for all lattices considered, regular 
or random, for all values of the rewiring probability q. 

As argued in [§| , the existence of this resonance effect 
is the result of a degradation of order. In our case, the 
degradation of order has its origin in the increasing im- 
portance of the inhibitory connections. This is clearly 



seen in figure [3] where we plot the standard order pa- 
rameter m = (ni(i)) as a function of the probability p of 
inhibitory links. The optimal probability for resonance 
(location of the peak of figure [T|) is found near the phase 
transition between the ferro and paramagnetic regions. 




100 200 300 400 500 600 700 800 900 1000 1100 
t 



FIG. 2: Evolution of magnetization in time (random network, 
q = 1, k — 10). Other parameters are: A'^ — 100, a = 0.15, 

O — — K — 1 

" — 100' ^ ^■ 




FIG. 3: Modulus of the average magnetization as a function of 
the probability of repulsive links. Same parameter values and 
symbol meanings than in in figure[T] In the regular networks, 
the existence of metaestable states reveals itself in a smaller 
magnetization at p = 0. 

The existence of this order-disorder transition and its 
relation to the resonance effects are reproduced by a sim- 
ple mean-field theory that we develop in some detail in 
the next section. 



MEAN-FIELD APPROACH 

At each time step the magnetization m(<) may change 
due to the modification of a single variable fj,i. The fol- 
lowing relation holds exactly for the ensemble average 
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m{t): 



In the case m > 0, 



Nm{t + t)= Nm{t) + {fi,{t + t) - (4) 

where — . . . , /iAr(i)) denotes the particular 

reahzation of the fii variables and (. . . | . . .) denotes a 
conditional ensemble average. By identifying r = 
and rearranging we get: 



m{t + r) — m{t) 



= {fidt + T)-fi,{t)\{fl{t)}) - 

-m(i) + (Ai.(i + r)|{Ai(i)}) (5) 



We now identify the left hand side as the time derivative 
and use the dynamical rules given by Eq.(IT]) to write: 



dm{t) 
dt 



= -m{t) + \f{t)\{sign[fm{^,{t)}) 



Design 



E 



where we have used the notation f{t) ~ asin(r2i). 
Since the forcing ](t) is independent of the state 
{/^}, then (sign[/(t)]|{^(t)}) = sign[/(t)]. Moreover 
|/(t)|sign[/(t)] = ]{t). For the last term of the right 
hand side of this equation we use the mean-field approx- 
imation: 



■ m(t) 



(7) 



where we replace the value by the average value 

■m{t). 

Now I]jey(i) "^y = K ^ '^K =^4^(1 + '^)^ a-i^d 
the mean-field approximation can be rewritten as: 



sign 



(-1) • Prob{[k+{l + k) - fc,K] m{t) < O) + 
(+1) • Fro5 ([fc+(l + k) - fc,K] m(<) > O) 

1 - 2Prob{[k+{l + hK] m{t) < O) 
= G{rn{t)) (8) 

from where we obtain the desired mean-field equation: 
^ = -mit) + /(<) + (1 - |/(<)|)G(™(t)) (9) 

The function G{m) can be easily computed in terms of 
the cumulative probability function Fk- of the binomial 
distribution of the number of positive links, given that 
the total number of links is fc^. This is precisely defined 
as: 



k + 



(10) 



Prob{ [k+{l + k)- k,K] m{t) < O) = 

„ J / 1 + kjK \ „ I kjK 
Prob kj < 



1 



1 



(11) 



while, for m < 0, 



Prob{[k+{l + k) - kiK] m{t) < O) = 



Prob kt > 



1 + K 



= l-Fi 



ki 



1 + K 



(12) 



By averaging over the distribution of the number of 
neighbors, we get: 



G(m) = sign(m) ^ Pk 



l-2Ffc 



1 + K 



(13) 



Pfc; being the probability that a node has ki links. Within 
the spirit of the mean-field approximation we assume that 
all nodes have the same number of links ki = k and 
replace the above formula by: 



G{m) = sign(m) 



k-K 



(14) 



k+<x 



In case of no forcing, f(t) = 0, the equilibrium value 
Too of the magnetization satisfies /tiq = G(mo). A stan- 
dard analysis of this equation predicts a phase transition 
separating a regime of non-zero stable solutions ±too ^ 
from a regime in which the only solution is toq = 0. The 
coexistence line is mo = 1 — 2Fk (^i^) and the critical 

point occurs at Fk (^i^^ ~ 1/2. In figure 0] we plot 
the equilibrium magnetization mo as a function of the 
probability p for fixed k. It is clear from this figure that 
the mean- field approximation reproduces the loss of order 
that arises as the proportion p of negative links increases, 
although the precise location of the transition point is not 
well reproduced. 

In figure O we plot the amplification factor computed 
after a numerical integration of Eq.®. Qualitatively, 
the results agree with those of simulations presented in 
the previous section: there is a resonance effect, i.e. the 
response shows a maximum as a function of p. The max- 
imum value is reached for a value Pc, close to that signal- 
ing the order-disorder transition. Furthermore, it can be 
noticed that the size of the amplification region, defined 
as the set of values of p for which i? > 1, is similar to 
the size of the transition region, defined roughly as the 
set of values of p for which the magnetization satisfies 
m{p) < 0.5 and the maximum is achieved at a value of 
p such that m{p) « 0.2 — 0.3. As the average number of 
neighbors k increases, the size Ap of this region decreases 
as and it disappears in the limit k — > oo. Since 

the relative dispersion in the number of positive links 
also scales as cr[fc+]/(fc+) fc~^/^, one is tempted to at- 
tribute the existence of the resonance to the existence of 
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FIG. 4: Modulus of the average magnetization as a function of 
the probabiUty of repulsive links according to the mean-field 
theory for k = 1. 

such a dispersion, a fact already stressed in the study of 
synchronized oscillations induced by diversity This 
is supported by a modified version of the mean-field ap- 
proach in which the dispersion is strictly equal to 0. This 
can be achieved by using in (|14|) the probability distribu- 
tion that would arise if all nodes had the same number 
of positive links, namely Fk{x) — Q ii x < pk and 
Fk{x) = 1 if X > pk. As shown in figure O in this case 
the amplification region has disappeared altogether. 
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FIG. 5: Spectral amplification factor versus probability of 
repulsive links according to the mean-field theory for a = 0.15, 

O — — K — 1 

" — 100' ^ 



MECHANISM 
Microscopic point of view 

We now give an explanation of some features of the 
observed resonance from a microscopic point of view, i.e. 
analyzing the evolution of individual values of fii. 



According to the rules ([T]), a chosen node takes the 
sign of the external signal with a probability |asin(r2t)|, 
independently of the current system configuration. To 
enhance resonance, there are two necessary requirements 
after a node has changed its state: to maintain the per- 
turbation in the next time steps, and to spread it to its 
neighbors. The crucial issue is then how the local config- 
uration of nodes and links helps (or hinders) this ordering 
process. 

To spread a perturbation, it would be an advantage 
to have all-attractive couplings; however, to maintain its 
state, the node cannot be too constrained by its neigh- 
bors. With a high homogeneity of the neighbors states 
and a positive connection with all of them, a perturbed 
spin would likely be forced to go back to its original state 
next time it is selected. At the other extreme, when all its 
connections are negative, a perturbed node is also very 
much constrained by the state of its neighbors, the local 
field being maximal for a local anti-ferromagnetic order- 
ing. At an intermediate level of positive and repulsive 
connections, we have the optimal state. It has a capac- 
ity to spread a perturbation to the whole network, but 
constrains minimally a node that has been perturbed. 
Due to the combination of attractive and repulsive links, 
the local field around a node is close to zero. Therefore, 
if a node changes its state, it possibly won't be forced 
to return to its previous position after consulting with 
its neighbors. On the other hand, it is easy to spread a 
perturbation: if a node had previously a zero local field, 
after one neighbor has changed, the balance is broken, 
and it has to align with that neighbor, if the connection 
is positive. 

This microscopic picture will help us to understand 
some of the observed features. For example, in figure [6] 
we show that the amplification region Ap decreases with 
increasing k whereas Pc tends to 0.5. Both facts agree 
qualitatively with the predictions of the mean-field the- 
ory. It is clear that for large k the condition of a local 
field close to zero can only be satisfied for a probability 
of repulsive links near 0.5. This is easily illustrated when 
one considers the case of p far from 0.5 and a uniform 
magnetization (at the peak of a signal's cycle). Getting 
a local field close to zero when the connectivity is high 
requires many neighbors flips. Since the unit to be up- 
dated is chosen randomly at each time step, it is likely 
that a unit is chosen twice before enough of its neigh- 
bors have been perturbed. On the other hand, p — 0.5 
is the upper limit for the amplification region, because 
a majority of positive links is necessary to have pertur- 
bation spreading. As the proportion of repulsive links 
approaches 0.5, more neighbors have a negative connec- 
tion and they will exert, when perturbed, an influence 
opposite to the signal. 

Note that for the resonance to disappear we need for- 
mally the limit k oo. In a finite network, the maximum 
value is A: = — 1 and, as shown in figure [5] for A^ = 201 
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and N = 1001, the resonance does not disappear com- 
pletely even for the maximum connectivity. 




P 



FIG. 6: Spectral amplification factor versus probability of 
repulsive links for a random network with q = 1, a = 0.15, 
— j^, K = 1. Main graph uses TV — 1001 while the inset 
shows the case A'^ — 201. 
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FIG. 7: Spectral amplification factor versus probability of 
repulsive links for a "no-dispersion" network in which all sites 
have the same number of positive and negative links. Due to 
the particular way the network is constructed [2^, only values 
of p = k/N where the total number of neighgbors per site, 
k, is an even integer number are allowed. Parameters are 
iV = 1001, a = 0.15, n = K = 1 (main graph) and 

= 201 (inset). Note that the amplification region shrinks 
as k increases. For comparison, we also include as dotted lines 
the results of figure [6l 



As we did in the mean- field treatment, and in order 
to isolate the influence of competitive interactions from 
the disorder induced by the dispersion in the number of 
links, we also present in figure [7] results from random 
networks when all nodes have exactly the same number 
of neighbors k and the same proportion p of repulsive 
links [26|. At variance with the previous results, an al- 
most total reduction of the amplification region can be 
achieved even for finite values of iV, for large enough k. 
This proves that diversity in the number of positive links 
is an important ingredient for the resonance effect. 

Why does dispersion matter? The precise mechanism 
is hard to grasp, but it is certainly related to a degra- 
dation of order at local level. To decrease the chance of 
having perturbed neighbors driving several units in the 
direction opposite to the signal, there have to be many 
nodes with a clear majority of positive links. But as we 
saw above - assuming every node had the same number 
of negative links - those units require many neighbors 
flips, to maintain their local field close to zero. However, 
if the nodes are heterogeneous, an unit with a lower than 
average number of repulsive links can profit from those 
neighbors that have many negative connections to other 
nodes. Since those are more susceptible to changes, their 
presence decreases the local field, thereby diminishing the 
need for many neighbors updates. This result confirms 
the importance of diversity in making the phenomenon 
more robust, but also shows that we can have an ampli- 
fication even without diversity. 



Macroscopic point of view 

In this subsection, we consider the explanation of the 
resonance from the macroscopic point of view, i.e. we 
look at the behavior of the collective variable (magne- 
tization) m(t). We assume that the dynamics of this 
macroscopic variable in the no-forcing case, / = 0, can 
be described in terms of relaxation in a potential func- 
tion V{m). The absolute minima ±mo of the poten- 
tial give the rest states which are separated by a po- 
tential barrier AV . This picture has proved to be valid 
in other problems with diversity in the parameters 
and it certainly holds in the mean-field limit where, ac- 
cording to the previous section, the dynamical equation 
is IF = with a potential V{Tn) = ^ - M{p)\Tn\ 

with M{p) = 1 - 2Fk (iq^) running from M(0) = 1 
to M{1) = —1. There are two minima of the poten- 
tial, mo = ±M{p) for M{p) > 0, and a single min- 
imum mo = for M{p) < 0, or p > pc, the critical 
point. For small p the barrier separating the two minima 
AV ~ MML- is high and it can not be overcome by the 
effect of the weak forcing f{t). The only effect of the 
forcing is a small oscillation around one of the minima 
(chosen by the initial conditions). As p increases, the two 
minima of the potential get closer to each other and the 
barrier separating them decreases such that, at a partic- 
ular value of p the forcing is able to overcome the barrier 
and m{t) oscillates between the two minima ±mo. As 
p crosses the critical value Pc, the two minima merge at 
mo = 0, the barrier disappears and the effect of the fore- 
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ing is reduced again to small oscillations around a single 
minimum. 

To apply this potential image beyond the mean-field 
approximation we need to include an important modifi- 
cation. As discussed before, the energy landscape is that 
of a spin-glass with many metastable states and two ab- 
solute minima ±mo. As a consequence, in the no-forcing 
case, the final state reached depends strongly on initial 
conditions. This is illustrated in figure [8] where we plot 
the probability distribution of the final magnetization. 
If the initial state is the ordered state fii = +1 (resp. 
— 1) Vi, the final magnetization is peaked near m = +1 
(resp. —1). If the initial state adopts Hi = ±1 randomly, 
then the final magnetization is peaked around to = 0. 
This reflects the existence of many barriers separating 
the metastable states from the absolute minima of the 
potential. When the forcing is introduced, it has to be 
able to overcome all these intermediate barriers. The 
final image is that of a particle moving in a "rugged" 
potential. As p increases, the height of those barriers de- 
creases and the forcing is able to explore a larger fraction 
of the configuration space, but not necessarily leading to 
trajectories ending in the absolute minima of the poten- 
tial. This can be seen in figure [5] where we show the 
effect of a forcing weaker than that used in figure [H The 
magnetization oscillates around a mean value that drifts 
with time. If we enlarge the period of the forcing - fig- 
ure [10] - the oscillations become wider and the system 
has now enough time to reach the equilibrium minima 
close to Too = ±1. p?} 
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FIG. 9: Evolution of the magnetization following a weak sig- 
nal a — 0.08, fl — j^, in the case of a random network, q — 1, 
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FIG. 10: Same as figure [9] for a slower forcing Q — glf-g. 
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FIG. 8: Distribution of stable states at the optimal probabil- 
ity pc — 0.25 in the case of an unforced random network with 
q — 1, N = 100, k = 10, K = 1 starting from three different 
inital conditions: all spins equal to -fl (data set indicated as 
m{t = 0) = 1), all spins equal to —1 {m{t = 0) = — 1) and 
spins take randomly the value ±1 ({m{t = 0)) = 0). 



CONCLUSION 

We have used Monte Carlo simulations and analytical 
(mean field) calculations to investigate the response of a 
system of two-state units, with both attractive and repul- 
sive interactions and majority-rule dynamics, to a weak 
periodic signal. For both regular and random networks, 
we have found that competing interactions can enhance 
the system response ~ a kind of "divide and conquer" 
strategy. In each case, a resonance was found for an op- 
timal percentage of negative links which depends on the 
model parameters. Applications include opinion dynam- 
ics and neuron networks but the model is generic enough 
to predict that the same type of effect can be found in 
other systems. 

We have discussed in some detail the microscopic 
mechanism for the amplification. We argued that the 
flexibility of the system to follow the external signal re- 
quires that the local field seen by each unit is kept close to 
zero and analyzed how this condition might be achieved 
in some parameter limits. 
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A macroscopic analysis, in terms of a relaxation dy- 
namics in a bistable potential, is able to explain the 
mean-field results. It is difficult to use this descrip- 
tion beyond the mean field treatment, due to the pres- 
ence of many metastable configurations. Because of their 
presence, a large response, corresponding to oscillations 
around (symmetrical) absolute minima can be obtained 
for a sufficiently slow forcing. There are studies that 
point to the role network topology plays in synchroniza- 
tion or response to stimuli [28|. Analyzing the eff'ect of 
coupling strength, degree distribution and other network 
characteristics on the coherent response may shed some 
light on how the mechanism can be optimized. 
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